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I.  Introduction 

The  deformation  and  failure  of  fiber  reinforced  composites  are  controlled  by  a 
variety  of  damage  mechanisms  occurring  at  the  fiber/matrix  interface,  in  the  matrix  and 
fibers,  and  among  fibers  in  the  fiber  bundle.  The  goal  of  this  program  was  to  develop 
accurate  micromechanical  models  of  these  damage  mechanisms  to  provide  predictive 
models  of  macroscopic  composite  behavior.  Such  predictive  models  clearly  highlight 
the  key  controlling  mechanisms  of  various  stages  of  deformation  and  to  quantitatively 
relate  the  constitutive  properties  of  the  fibers,  matrix,  and  interface  to  the  overall 
deformation  of  the  composite.  The  understanding  gained  through  such  predictive 
models  is  then  useful  for  optimizing  material  design  and  properties  for  various 
composite  applications  and  for  providing  a  baseline  for  studying  issues  such  as 
composite  degradation  due  to  creep  and  time-dependent  damage.  Considerable 
progress  was  made  toward  accomplishing  these  goals  during  the  grant  period;  much  of 
the  work  has  been  published  in  the  literature,  and  references  are  provided  at  the  end  of 
this  report.  Of  particular  relevance  are  two  major  invited  publications  (References  1 
and  2)  that  summarize  much  of  the  work  of  the  Principal  Investigator  in  this  field. 
Here,  we  only  summarize  our  progress  in  understanding  three  different  phenomena  in 
fiber-reinforced  composites:  the  tendency  of  matrix  cracks  to  deflect  at  a  fiber/matrix; 
prediction  of  the  stress-strain  behavior  of  a  ceramic  composite;  and  the  dependence  of 
tensile  strength  on  fiber  load  transfer. 

II.  Fiber/matrix  Interface  Debonding 

Consider  a  matrix  crack  in  a  CMC  that  approaches  the  fiber/matrix  interface 
(Figure  1).  If  the  matrix  crack  penetrates  and  fails  the  fiber,  then  the  composite  is 
brittle.  If  the  crack  deflects  at  the  interface,  then  the  strength  of  the  fibers  can  be 
accessed  and  strong,  damage  tolerant  behavior  can  be  obtained.  Understanding  the 
factors  controlling  deflection  versus  penetration  is  thus  critical  to  the  material  design  of 
CMCs.  A  crack  will  advance  in  a  particular  mode  if  the  energy  release  rate  G  due  to 
the  crack  advance  is  sufficient  to  pay  the  energy  cost  T  of  creating  new  crack  area.  The 
tendencies  for  deflection  and  penetration  can  thus  be  assessed  by  determining  the 
respective  energy  release  rates  Gd  and  Gp  and  comparing  these  to  the  respective  fracture 
energies  Ti  and  Tf.  At  the  interface  between  two  materials  of  differing  moduli, 
however,  the  energy  release  rates  diverge  or  approach  zero  as  the  amount  of  crack 
advance  (ap  into  the  fiber  or  ad  along  the  interface)  approaches  zero.  Previous  work  by 
He  and  Hutchinson  (HH)  circumvented  this  issue  by  considering  the  ratio  of  Gd/Gp  with 
ap=ad.3  Their  results  for  Gd/Gp  versus  the  Dundur's  parameter  a  (essentially  (Ef- 
Em)/(Ef+Em))  are  shown  in  Figure  2;  deflection  is  predicted  for  Ti/Tf  <Gd/Gp.  However, 
for  finite  ad,ap  and  finite  fiber  volume  fraction  Vf,  the  Ga/Gp  can  be  rather  different  than 
obtained  by  HH. 

We  have  performed  a  set  of  calculations  for  the  energy  release  rates  Gd  and  Gp  at 
finite  Vf,  finite  ad,ap  and  in  an  axisymmetric  geometry.4  The  calculations  are  performed 
using  a  model  due  to  Pagano  et  al.,  which  is  an  accurate,  efficient  method  for 
determining  stresses  and  energies  for  axisymmetric  composites.5  For  small  Vf=l%  and 
small  ad=ap=0.002r  (r=fiber  radius)  our  results  for  Gd/Gp  are  essentially  identical  to 
those  obtained  by  HH.  However,  at  V f=40%  and  crack  sizes  ad=ap  ranging  from  0.002r 
to  0.025r,  we  find  much  lower  ratios  of  Gd/Gp,  as  shown  in  Figure  2.  Crack  deflection 
is  thus  more  difficult  (requires  less-tough  interfaces)  as  both  Vf  and  aj=ap  increase 


from  zero ;  the  HH  prediction  tends  to  be  too  optimistic  in  assessing  prospects  for  crack 
deflection. 

A  group  at  WPAFB  has  fabricated  model  composites  composed  of  SCS-0  SiC 
fibers  and  a  glass  matrix.6  From  the  established  matrix  and  fiber  moduli  and 
toughnesses,  both  a  and  the  ratio  of  Ti/Tf  can  be  determined,  as  shown  in  Figure  3.  The 
HH  criterion  predicts  crack  deflection;  our  calculations  predict  crack  penetration  for  a 
wide  range  of  possible  crack  sizes  a<i=ap.  The  physical  composites  show  brittle 
behavior  and  no  evidence  of  fiber/matrix  interface  debonding.  The  observed  behavior 
is  thus  consistent  with  our  calculations  (although  our  minimum  values  of  a<i,ap  are 
determined  by  computational  limitations  not  physical  lengths). 

Our  results  show  a  fairly  strong  dependence  on  the  crack  sizes  a<j  and  ap  and  raise 
some  important  conceptual  issues.  How  are  ad  and  ap  physically  interpreted?  Are  they 
atomistic  lengths,  lengths  of  some  microstructural  scales  in  the  interface  and  fiber,  or 
lengths  of  intrinsic  cracks  in  the  interface  and  fiber?  The  interpretation  is  uncertain  but 
has  an  important  implication  for  interface  design.  Specifically,  if  the  conditions  for 
interface  debonding  depend  on  some  pre-existing  physical  defects  then  the  engineering 
of  interfaces  must  focus  on  controlling  such  defects.  This  is  far  more  difficult  than 
controlling  the  intrinsic  fracture  energies  of  potential  interface  coatings.  More  work 
must  be  done  to  properly  resolve  these  issues. 

III.  Stress-Strain  Behavior  in  CMCs 

After  debonding  has  successfully  occurred  at  the  fiber/matrix  interface,  the 
subsequent  deformation  of  the  composite  is  controlled  by  a  residual  sliding  resistance  x 
across  the  interface.  Furthermore,  additional  pre-existing  matrix  flaws  can  grow  and 
become  full-length  matrix  cracks,  and  pre-existing  flaws  in  the  brittle  fibers  can  also 
extend  across  individual  fibers.  These  cracking  modes,  and  the  associated  slip  at  the 
fiber/matrix  interface,  lead  to  non-linear  stress-strain  behavior  and  ultimately  to 
composite  failure.  Since  the  macroscopic  c-e  behavior  is  used  for  engineering  design, 
it  is  important  to  fully  understand  its  dependence  on  the  properties  of  the  fibers,  matrix, 
and  interface. 

Both  the  matrix  and  fibers  are  intrinsically  brittle  materials,  so  their  failures  by 
cracking  are  determined  by  statistical  distributions  of  flaws.  The  number  of  flaws  per 
unit  length  combined  with  the  interface  sliding,  which  controls  the  amount  of  strain 
released  per  crack,  determines  the  detailed  shape  of  the  stress-strain  curve.  We  have 
recently  developed  models  to  predict  the  evolution  of  statistical  matrix  cracking, 
statistical  fiber  cracking  and  the  associated  stress-strain  behavior,  as  a  function  of  the 
fiber  and  matrix  strength  distributions,  the  elastic  moduli,  volume  fractions,  and 
interfacial  x  in  unidirectional  CMCs.7'9  The  most  recent  models  demonstrates  an 
important  transition  in  failure  behavior.8,9  For  matrix  strengths  larger  than  the  fiber 
bundle  strength,  there  are  few  matrix  cracks  prior  to  composite  failure  and  hence  very 
little  non-linear  deformation  and  a  small  failure  strain.  For  lower  matrix  strengths, 
more  extensive  matrix  cracking  occurs  prior  to  failure,  allowing  for  increasing  non¬ 
linear  deformation  and  larger  failure  strain.  The  failure  strain  increases  until,  for  the 
weakest  matrices,  the  failure  strain  is  controlled  entirely  by  the  fibers.  Predictions  of 
this  transition  with  decreasing  matrix  strain  are  shown  in  Figure  4  for  a  "typical" 
composite.  Many  of  the  micromechanical  parameters  entering  the  theory  depend  on  the 
interfacial  sliding  resistance  x.  The  model  provides  predictions  of  the  stress-strain 


behavior  expected  as  t  is  varied  while  other  parameters,  such  as  fiber  and  matrix 
statistical  flaw  distributions  and  elastic  constants,  are  held  fixed.  Changing  x  can  occur 
due  to  modified  interface  coatings  or  interface  roughness  with  no  other  changes  to  the 
composite,  and  so  is  a  feasible  means  of  composite  design.  Figure  4  shows  the 
predicted  stress-strain  behavior  for  a  "tough"  composite  as  x  is  increased  by  multiples 
of  2,  4,  and  8  over  a  baseline  value.  Although  the  tensile  strength  inceases 
monotonically  with  increasing  x,  the  composite  failure  strain  does  not  increase.  In  fact, 
the  failure  strain  can  be  non-monotonic.  Such  unusual  behavior  was  previously 
unanticipated.  From  a  design  perspective,  it  is  important  to  have  both  high  matrix 
cracking  stresses  and  large  failure  strain  in  CMCs.  The  present  model  shows  the  trade¬ 
offs  maximizing  these  features  of  the  deformation. 

Application  of  the  model  to  SiC/SiC  "minicomposites"  tested  by  Lissart  and 
Lamon10  shows  excellent  agreement  with  experiment  when  only  the  matrix  strength 
parameters  are  adjusted  (Figure  6).8  These  materials  have  fairly  high  matrix  strengths 
and  so  (i)  the  failure  strains  are  lower  than  the  maximum  possible  values  (which  is 
0.88%  in  this  case)  and  (ii)  matrix  cracking  is  not  fully  saturated  at  the  failure  point. 
From  the  matrix  strength  parameters,  fiber  pullout  lengths,  and  hysteretic  behavior, 
three  independent  values  for  x  can  be  derived  and  all  predict  x«20  MPa. 

The  broad  capabilities  of  the  model  for  predicting  stress-strain  behavior  have  yet 
to  be  fully  exploited.  We  anticipate  wide  application  of  these  results  to  various  CMC 
systems,  particularly  as  the  manufacturing  reliability  of  components  from  specific 
material  systems  becomes  an  issue.  In  such  applications,  the  model  may  aid  in 
identifying  the  source  of  relatively  tough  versus  brittle  behavior  as  a  function  of 
variations  in  manufacturing  processes. 

IV.  Ultimate  Tensile  Strength 

Most  models  for  tensile  strength  assume  that  the  load  lost  by  broken  fibers  is 
transferred  equally  to  all  other  fibers  in  the  sample  cross-section  (Global  Load  Sharing, 
GLS).  Although  these  models  successfully  predict  the  strength  of  many  ceramic 
composites  having  low-x  interfaces,  it  is  expected  that,  for  higher-x  interfaces  and 
MMC  and  PMC  materials,  the  load  transfer  is  much  more  local.  The  existence  of  Local 
Load  Sharing  (LLS)  to  any  degree  also  implies  that  composite  strength  is  sensitive  to 
locally-damaged  regions,  is  dependent  on  the  physical  composite  size,  and  is 
statistically  distributed.  These  features  make  composites  weaker  and  more  brittle,  and 
thus  must  be  clearly  understood  at  a  fundamental  level  to  accurately  predict  composite 
behavior  and  to  guide  the  design  of  composite  components. 

We  have  developed  a  new  numerical  model  for  calculating  the  tensile  strength  of 
a  fiber-reinforced  composite  under  LLS  with  the  spatial  extent  of  the  load  transfer 
being  an  adjustable  parameter. 11,12  Numerical  simulations  of  composite  failure  exhibit 
all  of  the  features  noted  above,  and  particularly  the  scaling  of  strength  with  composite 
size.11'17  From  our  simulation  results,  we  have  deduced  some  fascinating  underlying 
scaling  relations  that  actually  relate  the  LLS  problem  to  the  analytical ly-tractable  GLS 
problem.12  We  have  also  generalized  the  simulation  model  to  incorporate  "notches"  or 
pre-existing  cracks,  free  surfaces  and  edges,  stress  gradients,  and  "random" 
arrangements  of  the  fibers. 

Careful  analysis  of  the  LLS  results  combined  with  size-scaling  concepts  has  lead 
to  the  development  of  a  handy  analytic  model  for  predicting  size-dependent  composite 


tensile  strength  under  conditions  of  local  load  sharing.  Our  previous  work  has  shown 
that  the  characteristic  in-situ  strength  and  gauge  length  for  fibers  in  a  composite  are 
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where  a0  is  the  Weibull  strength  parameter  at  a  reference  gauge  length  L0  and  p  is  the 
"Weibull  modulus",  r  is  the  fiber  radius  and  t  is  the  interfacial  sliding.  The  characteristic 
length  Sc  is  the  statistical  analog  of  the  “ineffective”  length  that  has  been  widely  used  in 
many  theories  of  composite  strength.  All  lengths  and  strengths  in  the  theory  are  naturally 
normalized  by  these  two  values,  respectively.  Now  consider  a  composite  consisting  of  n f 
fibers  and  physical  length  L,  with  fibers  having  a  Weibull  modulus  p.  There  is  a  critical 
damage  cluster  responsible  for  failure  that  consists  of  ni  fibers,  where 


n,  -  4O3y0_1 28  ,  (2) 

and  having  longitudinal  length  0.45c.  The  strength  distribution  for  the  critical  cluster  is  a 
Gaussian  distribution,  and  the  mean  strength  pC  and  standard  deviation  y*‘  (normalized 

by  crc)  for  the  appropriate  sizes  «/  are  given  in  Table  I.  The  full  composite  can  then  be 
considered  as  a  collection  of  mn  “links”  of  the  basic  critical-size  cluster,  with 

m=L/0.4Sc  ;  n=n/ni  .  (3) 

The  cumulative  probability  of  failure  for  the  full  composite  is  predicted  to  be  a  Weibull 
distribution, 


H„m(.CT)  =  1-e 
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where  the  characteristic  composite  tensile  strength  aum  and  Weibull  modulus  p  are 
given  by 
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The  measured  composite  strength  in  non-dimensional  form  is 

~  foe®, +0-/)<Tm  (7) 


where /is  the  volume  fraction  of  fibers  and  crm  is  the  load  carried  by  the  matrix  at  failure. 
Figure  7  shows  a  comparison  of  the  strength  vs.  composite  size  as  predicted  by  Eqs.  2-7 


and  as  obtained  via  the  LLS  simulations;  the  agreement  is  excellent  over  many  orders  of 
magnitude  in  composite  size. 
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Table  I.  Critical  link  size  «/ ,  mean  link  strength  p]h  and  link  standard  deviation 

/**  ,  normalized  by  oc ,  for  various  fiber  Weibull  moduli  p.  These  parameters  are  used 
in  the  analytic  model  of  Eqs.  2-7. 

As  an  example  of  the  predictive  capability  of  the  model,  we  have  applied  the  LLS 
model  to  predict  the  tensile  strengths  of  various  Ti-MMCs.  The  physical  composites 
are  actually  sufficiently  small  that  direct  numerical  simulation  of  the  full  composite  is 
feasible.  Input  to  the  predictions  are  only  constitutive  information:  the  in-situ  fiber 
strength  distribution,  fiber  volume  fraction,  fiber  diameter  r,  interfacial  shear  stress  t, 
and  matrix  yield  strength  amy.  The  predicted  and  measured  strengths  for  Ti-1100 
reinforced  with  SCS-6  fibers  are  shown  in  Table  II,  where  the  data  are  from  the  work  of 
Gundel  and  Wawner.18  Excellent  agreement  is  obtained,  which  is  typical  of  results 
obtained  on  other  Ti-MMC  materials.1 5 
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Ramamurty  et  al.  have  recently  taken  the  analytic  theory  and  applied  it  to  the 
prediction  of  strength  in  Aluminum  matrix  materials  reinforced  Nextel  Alumina 
fibers.19  With  no  adjustable  parameters,  the  theory  predictions  are  shown  in  Figure  8. 
The  agreement  is  excellent,  both  in  absolute  magnitude  and  in  the  trend  of  decreasing 
strength  with  increasing  size.  Ramamurty  et  al.  also  demonstrated  that  the  strength 
fluctuations  at  fixed  composite  size  were  due  primarily  to  the  fluctuations  in  average 
fiber  volume  fraction  among  the  samples  tested.  The  size-scaling  of  strength  thus 
cannot  be  obtained  by  using  a  simple  Weibull  model  relating  the  strength  distribution  at 
a  single  size  to  the  average  strength  versus  size. 

We  have  also  recently  applied  the  model  to  predict  the  tensile  strengths  of  several 
graphite/epoxy  polymer  matrix  composites.1617  Proper  application  of  the  model 
requires  accurate  input  data  on  fiber  strengths  and  critical  gauge  lengths.  These  were 
obtained  for  the  PMC  materials  through  careful  analysis  of  single  fiber  composite  test 
data  available  in  the  literature  but  interpreted  using  a  previously-developed  model  of 
the  PI.  Table  III  shows  the  constituent  fiber  strength  properties  and  critical  length,  and 
the  predicted  composite  tensile  strength  and  failure  strain.  The  predictions  agree  with 
the  experiments  to  within  20%. 17  Part  of  the  discrepancy  may  be  due  to  the  use  of  a 
single  Weibull  model,  whereas  data  on  Carbon  fibers  suggests  a  more  complicated 
strength-length  scaling  relationship.  In  any  case,  size-scaling  of  the  composite  strength 
plays  an  important  role  in  obtaining  quantitative  accuracy:  using  a  Global  Load  Sharing 
model  would  yield  predictions  that  were  too  large  by  40%  or  more. 
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1.5 
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2225 
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5863 
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0.64 

1.6 
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0.0137 

0.0164 

AS-4/ 

J2 
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10.7 

124 

0.53 

1.3 

xlO7 

1830 

1998 

0.0147 

0.0161 
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J2 
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10.7 
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0.54 

2036 

0.0130 

0.0161 

T300/ 

Epincote 

5436 

■ 

232 

1.00 

1.2 

xlO5 

3479 

0.0124 

0.0149 

T300/ 

#3601 

5181 

378 

■ 

138 

0.59 

1871 

0.0114 

0.0136 

T300/ 

#3631 

5436 

136 

0.59 

1963 

0.0128 

0.0144 

Table  III.  Constitutive  and  composite  data  for  various  systems,  and  measured  and  predicted  tensile 
strengths.  Th(/)  denotes  prediction  using  physical  volume  fraction;  Th(/i.)  denotes  prediction  using 
effective  volume  fraction.  Data  in  parenthesis  for  T300  systems  is  obtained  using  estimated  thermal  stress 
corrections  in  the  s.f.c.  test.  Units  are  as  follows:  o;(MPa),  <5>.(pni),  £c(GPa),  cr, „v(MPa). 

To  summarize,  the  numerical  and  analytical  tools  we  have  developed  are  predictive  and 
demonstrate  the  detailed  nature  of  fiber  damage  evolution  expected  in  fiber-reinforced 
composites  depending  on  the  nature  of  the  load  transfer.  These  models  will  form  the 
basis  of  future  work  on  time-dependent  damage  evolution  and  notch  sensitivity. 


V.  Summary 

We  have  briefly  discussed  our  major  progress  in  understanding  deformation  and 
failure  modes  in  fiber-reinforced  composites.  The  models  developed  here  show  the 
subtle  and  specific  connections  between  micromechanical  constitutive  properties  and 
macroscopic  measured  quantities  such  as  tensile  strength  and  strain.  Assessment  of  the 
models  by  detailed  comparison  with  experiments  demonstrates  their  accuracy.  These 
models  thus  provide  guidelines  for  material  development  and  optimization  of  the 
fiber/matrix  interface  toughness,  the  fiber  and  matrix  flaw  sizes,  the  interfacial  sliding 
resistance,  and  the  load  transfer  among  fibers.  Furthermore,  the  models  of  quasistatic 
deformation  and  failure  form  the  basis  for  the  development  of  micromechanical  models 
for  time-dependent  evolution  of  damage  and  ultimate  failure  including  the  phenomena 
of  creep,  interface  evolution,  and  fiber  degradation  at  elevated  temperatures.  Our  future 
work  is  aimed  at  developing  models  for  accurate  life  prediction  of  fiber-reinforced 
composites. 
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Figure  2.  G<j/Gp  vs.  a;  deflection  occurs  if  H/Tf  >G<j/Gp 
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Figure  3.  Gd/Gp  vs.  a;  deflection  occurs  if  P/rf  >Gd/Gp.  Experimental  data  for  p/Tf 
shown  as  +. 


Strain  (%) 

Figure  4.  Stress-strain  behavior  versus  characteristic  matrix  strength.  Note  the 
decreasing  failure  strain  with  increasing  matrix  strength. 
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Figure  5.  Predicted  stress-strain  behavior  for  interfacial  t  varying  by  factors  of  2,  4, 
and  8  from  a  baseline  value.  All  constituent  fiber  and  matrix  properties  are  unchanged. 
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Figure  6.  Predicted  (symbols)  and  measured  (solid  lines)  stress-strain  behavior  for 
SiC/SiC  minicomposites. 
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Figure  7.  Median  composite  strength  versus  dimensionless  composite  size,  as  predicted 
by  the  analytic  model  and  as  measured  in  simulations,  for  fiber  Weibull  modulus  of  5. 


Figure  8.  Predicted  (solid  line)  and  measured  (points)  strength  of  AI/AI2O3  composites 
versus  dimensionless  size  (after  Ramamurty  et  al.). 


